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Abstract 



To understand the mechanism allowing for long-term storage of excess en- 
ergy in proteins, we study a Hamiltonian system consisting of several coupled 
pendula in partial contact with a heat bath. It is found that energy absorp- 
tion and storage are possible when the motion of each pendulum switches 
between oscillatory (vibrational) and rotational modes. The relevance of our 
mechanism to protein motors is discussed. 
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Proteins are among the most important biopolymers for living systems. Still, the ques- 
tion of how proteins work is yet unanswered. Recently a noteworthing experiment concern- 
ing protein motors was performed M. In this experiment, the working process of a single 
molecule was directly investigated. Its results suggest that proteins often store energy ob- 
tained from a reaction with ATP and use it for work later (e.g., for enzymatic reactions with 
other proteins). The interval for energy storage was found to sometimes be very long, up 
to the order of seconds, while typical time scales for normal vibrations are several picosec- 
onds. How can proteins store excess energy for such a long time, somehow overcoming the 
relaxation process toward thermal equilibrium? In order for a protein to store energy for a 
sufficient time, energy must be absorbed into a certain part of it, in accordance with its own 
dynamics. Furthermore, some characteristic type of dynamics is required to hold the excess 
energy without losing it to the surrounding aqueous solution. 

As a first approach to understand the working mechanism of proteins, we construct a 
Hamiltonian system (in partial contact with a heat bath), which absorbs and stores energy 
for some time span, in spite of its eventual relaxation to thermal equilibrium on a much 
longer time scale. In this Letter, we adopt a system consisting of several coupled pendula 
Pll, each of which possesses two modes of motion, oscillation and rotation. We clarify the 
characteristic dynamics necessary for energy storage in connection with the coexistence of 
these two modes, and also show that partial contact with the heat bath is necessary for this 
storage. The relevance of our results to protein motors is also discussed. 

A protein consists of a folded chain of amino acids that assumes a globular shape [|J. 
Main chains are accompanied by side chains arranged around them, with each side chain 
hanging on the main chain as a kind of pendulum. Some side chains are gathered within a 
globular shapes, and a certain assembly of them can play an important role for the function 
of the protein. As an abstract model for the angular motion of side chains in such a functional 
assembly, we choose a system of coupled pendula. In particular, we study the idealized case 
of N identical pendula equally coupled to each other. Here, the oscillation of the pendula 
corresponds to the vibration of the side chains. With this simple model as an example, 
we show that energy absorption and long-term storage are generally possible in a class of 
Hamiltonian systems. Our study is restricted to a classical mechanical description, since 
quantum mechanical effects are believed to be irrelevant to protein dynamics (except for the 
choice of the potential). 

Our Hamiltonian is given by 

N 2 N 

H = K + v = Y. p f+ E (i) 

i=l z i,j=l 

where pi is the momentum of the i-th pendulum. The potential V is constructed so that 
each pair of pendula interacts through their phase difference with an attractive force to align 
the phases [@,|J: 

V{9 U Oj) = ^y^i 1 ~ cos(2n(6 t - %))}. (2) 

Hence the evolution equations for the momentum p^ and the phase are given by 

1 N 

P,- = ^£MM*<-*i)), (3) 
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Bj = Vy 



(4) 



In this model, all pairs of pendula interact identically. If the protein were a straight chain 
without folding, it could be modeled by a one-dimensional chain of pendula. In reality, 
due to the globular shape, there is a global interaction between pendula. Although the 
assumed "mean-field coupling" with uniform strength represents an extreme simplification, 
the resulting model can capture some general characteristics of protein dynamics. 

For the (dynamical) function of proteins, it is essential that they exist in aqueous solution, 
which plays the role of a heat bath. Since the hydrophobic part of a protein molecule is 
segregated from the aqueous solution, contact with the heat bath is restricted. Accordingly, 
we define the equation of motion so that only some pendula in the system contact the heat 
bath: 

pi = ~df~ 1Pi + v 2 ^W (*^)- ( 5 ) 

Here the heat bath is described by the Langevin equation, in which T represents its temper- 
ature, 7 is a relaxation coefficient, and is a Gaussian random form satisfying = 
and (£i(£i)£j(£2)) — Sij5(t 2 — ti), with (•) as the temporal average. An important point here 
is that contact with the heat bath is restricted only to the few pendula satisfying i < i c . 
This restriction may be interpreted as allowing for interaction only of the hydrophilic part 
with the aqueous solution. Note that this restricted heat bath (i.e. partial contact with the 
heat bath) is sufficient for realizing thermal equilibrium. 

We here briefly review typical behavior of Eqs.(|3]) and (|J) in the conservative case 
Pj0J§'@I' w hhout coupling to the heat bath. When the total energy is small, all pendula 
vibrate almost synchronized (i.e. \6i — Bj\ and \pi — Pj\ are small for all i, j). The typical 
time scale of vibration is ~ 1 - 10, while there is a sort of collective motion, giving rise to a 
periodic change of the degree of synchronization, with a time scale of ~ 10 2 - 10 3 . As the 
total energy increases, rotational behavior of pendula begins to appear. For a sufficiently 
large total energy, almost all pendula exhibit rotational motion. Here, each pi changes slowly 
in time, with values of \pi~Pj \ relatively large. In this case the pendula rotate almost freely, 
since the correlation of their phase with those of other pendula cannot be maintained, and 
the interaction term cancels out even with a short time average, except in the rare situation 
that two momenta take very close values. The time scale of rotation is 0(1). 

In the medium energy regime, rotational motion of a single pendulum appears inter- 
mittently from an assembly of vibrating pendula. Once a pendulum starts to rotate, it 
typically continues to rotate over many cycles, which is longer than the typical time scale of 
each pendulum's vibration (and rotation). We emphasize here that the effective interaction 
for a rotating pendulum is much weaker than those for vibrating ones, as mentioned above. 

When the system size N is very large |||5|||, there appears a phase transition at the 
energy E c /N = 0.0190 (T c = 0.0127), where the vibrational mode is dominant in the 'solid'- 
like phase and the free rotation is dominant in the 'gas'-like phase. Hereafter, we refer to the 
temperature range around the transition point (0.01 ~ T ~ 0.03) as the medium temperature 
range. 

Now, we consider the Hamiltonian system in contact with the heat bath described by 
Eq.fl5|). First, we note that there are two time scales, 7 -1 and for the coupling with the 
heat bath, where T = ^ic/N gives the dissipation rate of the system. Below the time scale 
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T^ 1 , the dynamics are approximated by the Hamiltonian system, while the dynamics of the 
pendula corresponding to % < i c are governed the time scale 7 . On the other hand, for 
a sufficiently long time scale, the system arrives at thermal equilibrium, irrespective of the 
values of 7, i c > 1, and the temperature. Note that the canonical ensembles are the same 
for the i c <^i N and i c ~ iV cases. 

Let us discuss how the dynamics change with i c by fixing the value of T. When almost all 
pendula are in contact with the heat bath, the time scale to arrive at thermal equilibrium is 
determined completely by the value of T (— 7) , independently of the dynamical properties 
of the system. On the other hand, in the case of partial contact i c <C N, the time scale 
to reach thermal equilibrium depends on the inherent Hamiltonian dynamics. For instance, 
the fluctuations of the total energy relax more slowly at both low and high temperature. 
In the medium temperature case, chaos (measured by the maximal Lyapunov exponent) is 
stronger []5]||. This is expected to hasten the relaxation. 

In the medium temperature regime, there is a large fluctuation of the total energy. 
Indeed, the dynamics of the total energy are correlated with the distribution of energy to 
each pendulum. Here, deviation toward larger total energy is supported by the concentration 
of energy in one (or few) pendulum (see Fig.p |7]. 

As the next step, we discuss the energetic behavior of the system far from equilibrium. 
Consider a special enzymatic event such as ATP attachment or its reaction. With such an 
event, proteins are moved far from equilibrium to allow for the accompanying conformational 
change. We study how such a change is related to the absorption and storage of energy. 

Due to the conformational change undergone with such a 'reaction' event, some portions 
of the side chains are forced far from the previous equilibrium state. This situation is 
modeled by introducing an external potential to cause a change in the shape of some part of 
the system as a result of the reaction. Adopting the coupled pendulum model, we assume 
that only one pendulum (to be chosen as the iVth pendulum) interacts with this external 
potential. Here we choose this external potential as VN,ext{® — On), where O is a certain 
phase variable representing the external part. For the sake of simplicity, we assume that 
(O — N ) can be approximately written as O + P^n with the constant (fi determining the 
stable relation between the two parts and the time scale parameter f3. Choosing a cosine 
potential as a simple example, we add the potential term 

V N , ext = V {1 - cos(27r(^ JV + 0o))}, (6) 

where Vq is a constant relating to the strength of the interaction. 

Then, we assume that the 'reaction' event occurs at t = to, thereby changing the relation 
between the pendulum On and the external coordinate O. This event can be accounted for 
by a change in the value of 0o that forces the pendulum far from equilibrium. For example, 
we consider the change from (f) Q = to (f>o = \- Then, after to, the equation of motion for 
the iVth pendulum is given by 

1 N 

Pn = tt^ E sin ( 27r (^ - M) ~ 27r/3y o cos(27r/30 JV ). (7) 

An example of time development for t > to is shown in Fig.^. Here, to is chosen so that 
On equals when the iVth pendulum is suddenly put in a high energy state parameterized 
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by Vq. In the computation we used a system with N = 10 and i c = 1; i.e. only the first 
pendulum is in contact with the heat bath. 

Through the interaction among the pendula, one (or few) pendulum absorbs a large 
portion of energy in FigfJ. This pendulum, after the relaxation of the Nth pendulum, 
continues to rotate in isolation with a large amount of energy. Then, it is affected very 
little by the other pendula over a long time interval. This fact allows for long-term energy 
storage. The reason that such a phenomenon exists is that the interaction experienced by 
the pendulum with such high energy is quite weak, as mentioned above. 

We define the lifetime for energy storage as the interval between the relaxation time of 
the Nth pendulum and that of the total energy (see Fig.Q). The latter is determined by the 
time tu at which the total kinetic energy K decreases to NT/2, i.e., the value at thermal 
equilibrium. On the other hand, the relaxation time of the iVtli pendulum is tentatively 
defined as the time £r at which p%/2 + Vpf >ext /N decreases to T. The results, however, do 
not depend on the specific choices of these relaxation times. 

The distribution of the lifetime is shown in Fig.|3](a), where 200 random sequences of 
are sampled. For comparison, we show the distribution for two cases with the same value of 
T: i c — 1 with 7 = 10 -2 and i c = N(= 10) with 7 = 10~ 3 . It is noted that the distribution 
of the lifetime is quite different for the two cases. For the case of partial contact, i c — 1, 
the typical lifetime of the energy storage is very long and reaches 10 7 . This contrasts with 
values of approximately 10 2 - 10 3 for the case with full contact, i c = N. This result suggests 
that partial contact with heat bath is necessary for long-term energy storage. 

Although the thermal equilibrium properties are common in the above two cases, a 
large difference appears when the system is placed far from equilibrium. The pendula 
corresponding to i < i c , interacting directly with the heat bath, cannot rotate freely over the 
I/7 time scale (in contrast with the Hamiltonian dynamics) and the dynamics are replaced 
by Brownian motion due to the heat bath. In the case of partial contact, a long duration of 
rotation is possible for pendula corresponding to i > i c only if the pendula possess sufficiently 
large energy to remain far from equilibrium. The pendula there become free from the thermal 
effect communicated through the distinctively weaker interaction with the other pendula, 
and nearly follow purely Hamiltonian dynamics. Although the prototype of this mechanism 
is observed as a large fluctuation around equilibrium in the medium temperature regime, it 
works clearly far from equilibrium. 

The temperature dependence of the mean lifetime of energy storage is displayed in 
Fig.|3](b). For i c = N, the mean lifetime is found to be almost the same for every tem- 
perature. Here, the distribution of lifetimes itself depends only on the temperature. The 
relaxation process is determined by the heat bath as in the case of thermal equilibrium, in- 
dependently of the temperature. On the other hand, for i c = 1, the relaxation time depends 
strongly on the detailed dynamical properties of the system at each temperature. 

The mean lifetime of energy storage has a maximum around T = 0.01 - 0.02. In this 
temperature range, the dynamics are most strongly chaotic, leading to rapid absorption of 
energy, and to the coexistence of rotational and vibrational modes, which is essential for 
energy storage. Also, the mixing properties of chaos are expected to determine the ability 
to transfer the excess energy localized at the iVth pendulum. A larger absorption of energy 
leads to a longer relaxation time. 

Energy absorption and storage in the case of partial contact is observed for a sufficiently 
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large V , small 7, and within a moderate range of (3 (1CT 3 ~ (3 ~ 1CT 2 ). V must be large in 
order to supply the system enough energy to kick it far from equilibrium. Unless 7 is small, 
the characteristic features of the Hamiltonian dynamics are largely destroyed by contact 
with the heat bath, and the ability of the system to absorb energy from the iVth pendulum 
is greatly weakened. If (5 is out of the above range, the time scale characterizing the motion 
of the iVth pendulum is much different from that of other pendula, and for this reason the 
other pendula effectively cannot interact with the iVth pendulum and they cannot absorb 
its energy. 

What type of Hamiltonian is required for long-term energy storage? Consider a system 
near the region of a phase transition. Here, some elements possessing higher energy "melt" 
and come to have a distinctively weaker interaction with the elements in the 'solid'-like 
phase. Such differentiation is indispensable to the storage of a large amount of energy. 
Partial contact with the heat bath is essential to maintain the differentiation of states needed 
to store the absorbed energy. Our coupled pendulum model provides the simplest example 
for this behavior. 

The dynamic mechanism for energy storage presented here is simple enough to be realized 
in real protein motors. It is only necessary that the protein dynamics exists near a phase 
transition region with 'solid'-like (vibrational) modes and 'gas'-like (rotational) modes and 
that the folded structure prevents some part of protein from experiencing the random effect 
of the heat bath. These conditions are expected to be satisfied for many kinds of proteins 
with a sufficiently large size, in addition to the protein motors. 

The authors are grateful to T. Yanagida, Y. Ishii and members of the Single Molecule 
Processes Project, JST for many meaningful suggestions. They would also like to thank T. 
Yomo, T. S. Komatsu, S. Sasa, and T. Shibata for stimulating discussions. This research was 
partially supported by Grants- in- Aids for Scientific Research from the Ministry of Education, 
Science and Culture of Japan. 
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FIGURES 



FIG. 1. Large deviation of the total energy E = K + V at thermal equilibrium, accompanied 
by the occasional appearance of the rotational mode (See the region 50000 <time< 60000). This 
behavior is characteristic of the medium temperature. For i c <C N, the rotational mode is mainly 
taken by pendula without direct contact to the heat bath. T = 0.02, N = 10, i c = 1 and 7 = 10~ 2 . 

FIG. 2. Process of energy storage at T = 0.02 in the case of partial contact with the heat bath 
(i c = 1 with N = 10). The values of the parameters are 7 = 10~ 2 , (3 = 10 -3 and Vq = 2. (a) Time 
series of the momenta p^. The iVth pendulum, kicked at t = 0, reaches thermal equilibrium at tm, 
as indicted by A. The excess energy is transferred to another pendulum {i 7^ 1,N), as shown by 
B. (b) The time series of the total energy E. The inset is a close-up of the graph for time> 10000. 

FIG. 3. Lifetime for energy storage when N = 10, T = 10 -3 , (5 = 10~ 2 and Vq = 4. (a) 
Distributions of the lifetime log 10 (i_R — £ro) fo r energy storage at T = 0.02 (medium temperature). 
Each distribution is obtained from 200 samples. The unshaded distribution is that for the case 
of partial contact (i c = 1 with 7 = 10~ 2 ), while the shaded distribution is that for full contact 
(i c = 10 with 7 = 10~ 3 ). (b) Dependence of the mean lifetime (tji — too) for energy storage on the 
temperature T. The values are averaged through 200 samples. The solid line denotes the case of 
partial contact (i c = 1 with 7 = 10 -2 ), while the dotted line denotes that for full contact (i c = 10 
with 7 = 10~ 3 ). 
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N.Nakagawa and K.Kaneko Figure 1 . 




N.Nakagawa and K.Kaneko Figure 2. 




N.Nakagawa and K.Kaneko Figure 3. 



